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Abstract 


Quantifying the “capability” of a manufacturing process is an important initial step in 
any quality improvement program. Capability is usually defined in dictionaries as “the 
ability to carry out a task, to achieve an objective”. Process capability indices(PCIs) is de¬ 
fined as a combination of materials, methods, equipments and people engaged in producing 
a measurable output. PCIs which establish the relationships between the actual process 
performance and the manufacturing specifications, have been a focus of research in quality 
assurance and process capability analysis. Capability indices that qualify process potential 
and process performance are practical tools for successful quality improvement activities 
and quality program implementation. As a matter of fact, all processes have inherent sta¬ 
tistical variability, which can be identified, evaluated and reduced by statistical methods. 
Generalized Process Capability Index, defined as the ratio of proportion of specification 
conformance (or, process yield) to proportion of desired (or, natural) conformance. We re¬ 
view the process capability indices in case of normal, non-normal, discrete and multivariate 
process distributions and discuss the inferential aspects of some of these process capabil¬ 
ity indices. Relations among the process capability indices have also been illustrated with 
examples. Finally we also consider the process capability indices using conditional order¬ 
ing and transforming multivariate data to univariate one using the concept of structural 
function. 
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1 Introduction 


Process capability is an important concept for industrial managers to understand. The chal¬ 
lenge in today’s competitive markets is to be on the leading edge of producing high quality 
products at minimum costs. This cannot be done without a systematic approach and this ap¬ 
proach is contained within what has been called “Statistical quality control”. The segment of 
statistical quality control(SQC) discussed here is the process capability study. Process capabil¬ 
ity indices (PCIs) aim to quantify the capability of a process of quality characteristic to meet 
some specifications that are related to a measurable characteristic of its produced items. These 
specifications are determined through the lower specification limit (L), the upper specification 
limit (U) and the target value (T). A variety of such indices have been developed in last two 
decades. Following the work of Kane (1986), a number of articles have appeared introducing 
new indices and /or studying the properties of existing ones. Kaminsky et al. (1998) have given 
their critical comments on uses of these indices, and suggested a future measurement. Uses of 
process capability indices in administration and supplier certification have been discussed in 
Latzko (1985) and Schneider et al. (1995). Excellent reviews on them are given by Rodriguez 
(1992), Kotz and Johnson (1993, 2002), and Kotz and Lovelace (1998). In addition. Spiring 
et al. (2003) provide an extensive bibliography on process capability indices. It appears to 
be a general acceptance of the idea that PCIs can be used only after it has been established 
that a process is in “statistical control” (for example, by the use of control charts). This is 
reasonable if it simply requires that there be no irregular changes in quality level. However, 
there seems to be, in some quarters, an assumption that the measured characteristic should 
have a normal distribution (at least approximately), although it is difficult to see why a good 
industrial process must result in a normal distribution for every measured characteristic. The 
majority of the process capability indices discussed in the literature are associated only with 
processes that can be described through some continuous distributions of the characteristics 
and, in particular, normally distributed characteristics. The most widely used such indices are 
Cp [by Juran (1974)], Cpk [by Kane (1986)], Cpm [by Hsiang and Taguchi (1985)] and Cpmk [by 
Choi and Owen (1990), Pearn et al. (1992)] and their generalizations for non-normal processes 
suggested by Clements (1989), Constable and Hobbs (1992), Pearn and Kotz (1994-1995), Pearn 
and Chen (1995) and Mukherjee (1995). A number of new approaches to process capability 
analysis have been attempted by Carr (1991) and Flaig (1996). Another index that is given by 
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Boyles (1994) is also worth mentioning. Often, however, one is faced with processes described 
by a characteristic whose values are discrete. Therefore, in such cases none of these indices can 
be used. The indices suggested so far whose assessment is meaningful regardless of whether the 
studied process is discrete or continuous are those suggested by Yeh and Bhattacharya (1998), 
Borges and Ho (2001), Perakis and Xekalaki (2002, 2005) and Haiti et al. (2010). 

In this article, our objective is to look into the overview of successive development of the process 
capability indices from the very beginning and also the inferential aspects of these capabilities, 
not only in case of normal process distribution, but also for non-normal and in some cases for 
discrete process distributions. The article has been organized as follows. Section 2 deals with 
some notations, where as Section 3 consider the background of process capability indices which 
have been developed in the earlier stage. Section 4, 5 and 6 consider the same in case of non¬ 
normal process distribution, any general process distribution, whether the process is normal, 
non-normal or discrete process distribution and multivariate process distributions respectively. 
Section 6 is devoted to gives concluding remarks. 


2 Notations 


L 

U 

fix) 

F{x) 

F{x) 

UDL 

LDL 

Oil 

0.2 

a 

p 

Po 


Lower specification limit. 

Upper specification limit. 

Probability density (or mass) function of quality characteristic X. 
Cumulative distribution function of quality characteristic X. 

1 — F{x). 

Upper desired limit. 

Lower desired limit. 

P{X < LDL). 

P{X > UDL). 

Mean of the distribution. 

Median of the distribution. 

Standard deviation of the distribution. 

Process yield, Le, f{x)dx=F{U) — F{L). 

Desirable yield, F{UDL) — F{LDL) = 1 — a with a = ai + a 2 . 
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3 Background 


Process Capability analysis is a technique that has application in many segments of the product 
cycle, including product and process design, vendor sourcing, production or manufacturing 
planning and manufacturing. Quality of the products has always been a major concern for 
both consumers and producers. There appears to be a general acceptance of the idea that PCIs 
can be used only after it has been established that a process is in “statistical control”. This is 
reasonable, if it simply required that there be no irregular changes in quality level. However, 
there seems to be, in some quarters, an assumption that the measured characteristic should 
have a normal distribution, although it is difficult to see why a good industrial process must 
result in a normal distribution for every measured characteristic. The first PCI was developed 
by Juran(1974) and is defined as 

C, = (3,1) 

_d_ 

where, d = {U — L)/2. Note that Cp does not depend on process mean. A Cp value of 1 means 
that 99.73% of all individual items will be within specification. However, if there is only a 
slight change in the process mean or a slight increase in the process variation, a larger portion 
of items will be out of specification. Currently, Cp values of 1.33 or greater are used widely 
as industry standards. Some companies now require Cp values of 1.66 from their suppliers. 
A Cp of 1.33 indicates that 99.73% of products production is within (100%/1.33) = 75% of 
the specification, whereas a Cp of 1.66 shows that almost all product measurements fall within 
(100%/1.66) = 60% of the specification. For example, if the width of the specification is about 
20 and the standard deviation is 3, 


Cp = 


30 - 10 


18 


= 1.1 


indicating that not only can the process produce almost all production within specifications but 
also 100%/1.11 = 90% of the specification will contain essentially all the production measure¬ 
ments. Whereas Cp values represent the setting in which the process is centered between the 
specification limits, process are frequently operated with the mean other than halfway between 
the upper and lower specification limits. 

When the process mean is off-centered of the specification, the result is that one specification 
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limit (the closer to the process mean) becomes the focal point of the process capability calcu¬ 
lation. The modified capability index, represents the situation in which one specification limit 
is more likely to be exceeded. The Cpk is then introduced by Kane(1986) to reflect the impact 
of /i (process mean) on the process capability indices. 

Cpk = min ^ ^ ^ (3.2) 

[ 3cJ 3cJ 

3a J ’ 

where, M = {U + T)/2. The change in the denominator from six to three standard deviations 
is the result of the two one-sided quality concerns. For example, the Cpk for the situation 
described above is Cpk = mm[0.78,1.44] = 0.78. When Cp and Cpk are calculated for a process 
that is centered, the two calculation methods give the same number. Therefore, Cpk is generally 
preferred because it is not dependent on the process being centered. However, the disadvantage 
of using only the Cpk value is that it does not reveal whether a poor process capability is the 
result of a process having a large dispersion or an off-centered process. 

Therefore, both Cp and Cpk should be used to judge the process capability. If the index values 
are different, then this is a sign that the process mean is off-center and considerable improvement 
in the proportion of acceptable product can be made by centering the process mean. Cpk has 
been criticized because it does not measure the process deviation from target. A process may 
have a high Cpk even when the process mean is off-target and close to the specihcation limits 
as long as the process spread is small. 

As an alternative, process capability index, Cpm, given by Hsiang and Taguchi(1985), takes into 
account the influence of the departure of the process mean n from the target T, 



where /x is the process mean and T is the target value and E{.) denotes “expected value”. It 
is easy to verify that Cpm will possess necessary properties required for assessing process capa¬ 
bility. If the process variance( i.e., ) increases (decreases) the Cpm will decrease (increase). If 

the process drifts from its target value (i.e., if fj, moves away from T) then Cpm decreases. In 
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the case where both the process variance and the process mean change, the Cpm index reflects 
these changes as well. Usually, T = M ; if T 7 ^ M the situation is sometimes described as 
“asymmetric tolerances” (see Boyles (1994) and Vannman (19976,1998a)). The measure Cpm 
sometimes called the “Taguchi index”. There is also the hybrid index given by Pearn et al. 
(1992) 


a, 


pmk 



1 

y: 

min 


d- 

]y-M\ 


2,^E{^x-TY 


— L 

+ (^ - rf 


(3.4) 


Clearly, Cpmk is only meaningful when deviation from target is the main concern. This index 
is based on the quadratic loss function and, thus, should only be used when there is evidence 
of a quadratic monetary loss. There are an infinite number of possible loss functions, and in 
many practical cases, material outside specification limits may result in a total loss rather than 
a quadratic loss. Clearly Cp > Cpk > Cpmk and Cp > Cpm > Cpmk- The relation between Cpk 
and Cpm is less clear cut. From Equations (1) and (2) we have 
Boyles (1994) suggested index is 


Spk = 


-1 


2 \ a 2 V cj 


[Boyles, (1994)] (3.5) 


An enlightening view of relations among our basic PCIs can be obtained from studies of the 
superstructure PCIs” introduced by Vannman (1995). Let d = {U — L)/2 be the half length of 
the specification interval. Vannman(1995) defined a class of PCIs which is a generalization of 
the existing indices. It depends on two non-negative parameters u and v 


Cp{u,v) = 


d — u 


\V annman, 1995] 


(3.6) 


6\/ + v{y — T)2 

where y is the process mean, a is the process standard deviation, T is the target value and 
M = {U + L)/2 is the mid-point of the specification interval. By varying the values of u and v, 
one can easily see that (7^(0,0) = Cp (Juran (1974)), Cp{l,0) = Cpk (Kane (1986)), Cp{0, 1) = 
Cpm (Hsiang and Taguchi (1985)), and Cp{l,l) = Cpmk (Pearn, Kotz and Johnson(1992)). 
Moreover, when T = M and for u,v > 0, Cp{u,v) has the same interpretation as Cpm and 
Cpmk ■ 
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Spiring (1997) also defines a PCI, = Cp(0, w). However, in this definition, w is not 

necessarily a constant; it may be a function of In principle this allows \w{^ — T)/(t]^ in 

Vannman’s formula to be replaced by any function of — T)la. So, in effect, 


cl"’ = i+a((iflr)/,rp ■ [Spirins, 1997] 


(3.7) 


with a general choice of function g{.), though for practical purposes it should be a positive, 
increasing function. 


4 PCIs under the assumption of non-normal process distribu¬ 
tion 

As already noted, the “6” in Equation (3.1) has been associated with the idea that a normal 
distribution for X provides a satisfactory approximation. Of course, both practitioners and 
theoreticians realized that this would not always be the case, and some (at least) of the second 
group energetically busied themselves with the task of coming up with relevant information and 
suggestions. Some practitioners, on the other hand, have claimed that Cp need not be assessed 
on the grounds of direct relevance to properties of NC product, though it is not clear what 
other means of assessment are to be used. 

At a relatively early date, Clements (1989), in an influential paper, suggested that “6it” be 
replaced by the length of the interval between the upper and lower 0.135 percentage points of 
the distribution of X (this is for a normal N(|U, fx^)). The new PCI is 

Cp = --— [Clements, (1989)] (4-8) 

sl—a sa 

where is defined by P[A < taking a = 0.00135, so that ^i_a, Ca are the upper and 

lower 0.135 percentiles of the distribution of X. [For a N(//, cj^) distribution = ^+3 cj, = 

H — 3cj.] 

Clements (1989) suggested fitting a Pearson system distribution for X, in order to obtain the 
required values. Applications of this kind of method, with various assumed distributional 
forms, have been quite numerous since 1992. References include: Rodriguez (1992), Bittanti et 
al. (1998), Lovera et al. (1997) -all Pearson system; Castagliola (1996) -Burr distributions; Far- 
num (1996), Polansky et al. (1998), Pyzdek (1992) -all Johnson system; Padgett and Sengupta 
(1996) -Weibull and log-normal; Mukherjee and Singh (1997) -Weibull; Sarkar and Pal (1997) 
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-extreme value; Somervile and Montgomery (1996/7) -t, gamma, and log-normal; Sundaraiyar 
(1996) -inverse Gaussian. As mentioned above, Polansky (1998, 2001) uses a general kernel 
fitting method. 

An Index defined on Inter quantile range is given by 

' = [Mutteriee, (1995)1 (4.9) 

Yeh and Bhattacharya (1998) proposed use of a PCI based on the ratios of expected proportion 
NC to actual observed or estimated proportion NC. Let X be the random variable associated 
with the characteristic of a process under study. Here we assume that X is positive and has a 
continuous distribution. Let po be the desired proportion of non-conforming output and p be 
the actual proportion of non-conforming. Then 

a) Cp = — ; [Yeh and Bhattacharya, 1998] (4-10) 

P 

Yeh and Chen (1999) have extended this to multivariate cases. Another PCI, suggested by Yeh 
and Bhattacharya (1998), distinguishes between NC items for which X is less than L, and those 
for which X is greater than U. Let Oq and o-q be the expected proportions of non-conforming 
products the manufacturer can tolerate on the lower and upper specification limits, respectively. 
Let CiL = P{X < L) and ajj = P{X > U) be the actual proportions of non-conformance of the 
process. The question of whether the process is “capable” can then be answered by comparing 
aQ to ai and to ajj. Thus 

fan 

b) Cf = min — , —^ ; [Yeh and Bhattacharya, 1998] (4.11) 

[aL auj 

Here the index is directly linked to the probability of non-conformance. The computation 
involved in estimating this index is more intensive than the conventional PCIs. This perhaps is 
not a serious constraints in today’s computer technology. A new measure of process capability 
which is in one-to-one correspondence with the process non-conforming fraction tt. Namely 

C = [Borges and Ho, (2001)] (4.12) 

where vr = 1 — P(L < X < U) is the in-control process fraction defective. The index C has some 
interesting properties, (i) It is one-to-one correspondence with the in-control process fraction 
defective tt. Processes with the same fraction defective will be equally capable. Moreover, the 
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process capability will respond only to changes in the process fraction defective and not to 
changes in the distribution of the observed quality characteristic. 

(ii) If, under process stability, X is normally distributed with mean jj, = {L + U)/2 and standard 
deviation a, then 



= 1 - $ 


U-fi 


a 


and 

i.e., process capability is measured by the value of a centered normal process with the same 
in-control process fraction defective tt. This property makes C values easy to interpret and 
understand, (iii) It can be applied to discrete and continuous, uni- and multivariate quality 
indicators and no assumption on the distribution of the observed quality characteristic is made. 


5 PCIs under the assumption of any general process distribu¬ 
tion 


Now a PCI is introduced which overcomes the drawbacks of the standard indices discussed 
earlier. It is defined as the ratio 




pc 


1 -PO . 

1 — p ’ 


[Perakis and Xekalaki, 2002] 


(5.13) 


Where pq denotes the minimum allowable proportion of conformance, p is the process yield 
lying between L and U. The value of pq must be intuitively close to unity and depends on the 
nature of the examined process and the requirements of the customers. This index is used for 
process with unilateral as well as bilateral tolerances and can take into account the minimum 
acceptable process yield(if such quantity has been set). 

A generalized process capability index given by Haiti et al. (2010), which is the ratio of pro¬ 
portion of specification conformance (or, process yield) to proportion of desired (or, natural) 
conformance. Almost all previously defined process capabilities are directly or indirectly asso¬ 
ciated with this generalized one. 


F{U)-F{L) 
F{UDL) - F{LDL) 

p_ 

Po' 


(5.14) 
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The following theorem gives different forms of Cpy when X has uniform distribution. 
Theorem 5.1 If X follows uniform distribution over (a,b),a < b, then Cpy reduces to 


(a) Cp = when LDL = fj, — 3a and UDL = ^ + 3a [Juran (1974)]; 

(b) Cp = LDL = and UDL = Ci-f with k = P{X < ^if) [Clements 

(1989)]; 

(c) I = when LDL = and UDL = F~^{1 — 02 ) [Mukherjee 

(1995)]. 


When proportions of desired conformances are ai for lower tail and 1 — 02 for upper tail, 

then Cpy = , which for a normal process with LDL = n — 3a and UDL = // + 3cj is 

r — p 
^py ~ 0.9973 • 

Theorem O motivates us to formulate the generalized index Cpy given in (I5.14h . Let us 
examine the behaviour of the proposed index Cpyy for different values of p. Obviously, if p = po 
(for normal distribution, p = 0.9973), then Cpy = 1. If the process yield p is greater than pQ, 
then Cpy > 1 whereas, if p < po; then Cpy < 1 , and the value of the index approaches zero as p 
tends to zero. Thus, the smallest possible value of Cpy is zero. 

Now the situation is of importance when the process is off-centered, i.e, F{L) -\- F{U) 7 ^ 1, 
but the proportion of desired conformance is achieved. In that case they defined the index as 
follows: 


c 


pyk 


= mm 


|T(t/)-F(/ie) F{p,)-F{L) 


2 - «2 


2 *^1 


— min { Cpyu , Cpyl } 


(5.15) 


When Cpy and Cpyk are calculated for a centered process, they come out to be the same. 
Therefore, Cpyk is generally preferred because it is not dependent on the process being centered. 
If the index values are different, then this is a sign that the process median (mean in case of 
symmetric distribution) is off-centered, and considerable improvement in the proportion of 
acceptable product can be made by centering the process median. 


When «! = 02 = f) then 


c 


pyk 


min 


FjU) - h k-nL) \ 

1 ( 1 _«)’ 1 ( 1 _«) j 


d-\l-F{M)\ 
^(1 - a) 
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where F(M) = 

For a normal process and under the assumption of uniformity as in Theorem 15.11 Cpyk reduces 
to 


a 


pk 


u - n 

mm < -, 

3c7 ’ 

d- |;U-M| 

3cj 


^^-L \ 
3cj / 


where d = ■ 

It generally happens that process target T is such that F{T) = . ■£ ^ F(L)+F(f/) ^ 

the situation may be described as “generalized asymmetric tolerances” [Boyles (1994) and 
Vannman (1997,1998) have described by the term “asymmetric tolerances” when T ^ M = 
Under this circumstance, they defined 




pTk 


= mm 


'f{U)-F{T) F(r)-F(L)| 


(5.16) 


I 2 ~ *^2 2 ~ 

At this point, it would be interesting to note that CpTk is equal to Cpyk when T = |Ue- 
This generalized process capability index overcomes many deficiencies of the PCIs which have 
been already discussed and more or less all the existing PCIs are directly or indirectly associated 
with this index. And it can be used for process with unilateral or bilateral tolerances. 


6 PCIs under the assumption of multivariate process distribu¬ 
tion 

A more precise title for this Section would be “PCIs for Use When X is Multivariate”. Many 
of the PCIs in this group are not, in fact, multivariate. May be they should be, but writers 
have opted for construction of univariate PCIs, based on the multivariate distributions of X. 
Nevertheless we will term them all “multivariate PCIs” (MPCIs). References with a title 
including the word “multivariate” or “bivariate” are: Beck and Ester (1998); Bernardo and 
Irony (1996); Boyles (19966); Chan et al. (1991); Davis et al. (1992); Hellmich and Wolf 
(1996); Hubele et al. (1991); Karl et al. (1994); Li and Lin (1996); Mukherjee and Singh 
(1994); Niverthi and Dey (2000); Shariari et al. (1995); Taam et al. (1993); Tang and Barnett 
(1998); Veevers (1995,1998,1999); Wang et al. (2000); Wierda (1992,1993,1994a, 19946); and 
Yeh and Chen (1999). Multivariate situations are also discussed in the following references. 
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that do not indicate, explicitly, in their titles that this is so: Chan et al. (19886); Wang and 
Chen (1998/9); and Wang and Hubele (1999,2001). 

The univariate specification interval (L < X < U) is now replaced by a specification region. 
This may just be constructed from separate specification intervals; one for each variable Xi in 
X. The specification region is then the hyperrectangle 

V 

n {Li <Xi< Ui). 

i=l 

However, more complex regions may be used, reflecting perceived relations among the variables 
in X. These are of the general form 


L < g{X) < U. 

Often, L is zero. Possibly for mathematical convenience, g{X) is often taken as a monotonic 
function of the joint probability density function of X. Thus if X is assumed to have a multi¬ 
variate normal distribution, one might take 

g{X) = {X - gyj:-\X - g) (6.17) 

and regard an item as NC if g{X) > U. In this way we obtain the ellipsoidal specification 
region 

{X - gy^-^{X - g) <U. (6.18) 


An analogue of Cp is 

Volume of{{X - gyj:-\X - g) < U} _ /t/V 
Volume of{(A: - gyV-^{X - g) < R} ~ \r) ' 

where 

P[{X - gyV-^{X - g) < R] = 1 - p. 


If the distribution of X is multivariate normal then {X — gyT,~^(X — g) has a distribution 
with V degrees of freedom and R = Xv,i-p upper 100(1 — p)% point of the x^ or “chi- 
squared” distribution with v degrees of freedom). 

Chen (1994) applies this method to the case when the specification region is of the form in 
Equation (I6.18|] . The region is defined by 


?7iCIXj=l^2,...,D 



< 1 
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with Mi = {Li + Ui)/2 and di = {Ui — Li)j2, and Chen defined MCp as R where 


P 


TnClXi=i^2,...,v 



< R 


1 — p. 


( 6 . 20 ) 


There can be many variants on these approaches. For example, the g{X) in Equation (I6.17P 
might be replaced by {X — p)'A~^{X — p) where ^4 is a positive matrix, not necessarily the 
variance-covariance matrix of the distribution of X. 

Shariari et al. (1995) proposed a truly multivariate MPCI. It contains three components. The 
first is of the type in Equation (I6.20p . The second is the significance level of the Hotelling’s 
statistic 


= n{X - pyS-\X - p), 


which is 


P 


Pv.n—v P 


n — V 


v{n — 1) 

where F^^n-v denotes a variable having the F distribution with v,n — v degrees of freedom. 
The final component just takes values 1 or 0 according to a modified process region -defined as 
the smallest region similar in shape to the specification region, circumscribed about a specified 
probability contour (of the distribution of X) -is or not entirely contained in the specification 
region. 

Wang et al. (2000) compared this 3-component MPCI with Chen’s (1994) MCp and with an 
index MCpm proposed by Taam et al. (1993) which is also a ratio of two volumes. The volume 
in the denominator is the same as in Equation (j6.19p with R = Xv,i-p with p = 0.0027, while 
in the numerator we have the volume of a “modified specification region” which is the largest 
ellipsoid centered at the target that is within the original specification region. 

Bairamov (2006) introduced the conditionally ordered order statistics for multivariate observa¬ 
tions. Let Xi, X 2 , ..., Xn G S C i?™ be i.i.d. random vectors with m-variate c.d.f. F{x) and 
p.d.f. f{x), where x = (xi, X 2 , ■■■, Xm) and S is the support of X. Consider the real-valued 
function N{X) : R^ R, x = (xi, X 2 , ■■■, Xm), which is continuous in its arguments satisfying 
N{X) > 0 for all x G T?™ with N{X) = 0 if and only if x = 0, where 0 = (0, 0, ..., 0). N{Xi), 
N{X 2 ), ..., N{Xn) are i.i.d. random variables with c.d.f. P{N{Xi) < N{x)), x € R. The 
function N{x) introduces partial ordering among the random vectors Xi, X 2 , ..., Xn and Xi is 
said to be conditionally less than X 2 (or Xi precedes X 2 ) if A^(Wi) < N{X 2 ). This ordering is 
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denoted by Xi -< X 2 . 

Selection o f the function N{x) is independent of the structure of the system and it can be cho¬ 
sen in various ways depending on the conditions available in the process/system. The following 
selections for N{x) may be of special interest. 

(a) N{xi, X 2 , Xn)= YA=lO‘i^i- 

(b) N{xi, X2, Xn)= min{xi, X2, ■■■, Xn)- 

(c) N{xi, X2, Xn)= max{xi, X2, x„). 

While the first function reflects the combined effects with the weights ai, the last two functions 
consider the extreme values. 

For multivariate quality characteristic, if there are lower and upper specification vectors L and 
U_ given, then L ^ X ^ is more meaningful and similarly for lower desirable vector LDL 
(or lower tolerance vector LTD and upper desirable vector UDL (or upper tolerance vector 
UTL ). In that circumstance Maiti et al. (2013) defined generalized process capability index 
for multivariate process quality characteristic as 


^ PjL^X^ U) 

Py P iLDL PX< UDL) 

P{N{L) < N{2Q < N{U)) 

PiN iLDL ) < N(X) < N (UDL )) 
F{N{U)) - F{N{L)) 

FiN iUDD ) - FiN iLDD ) 

= ^ 

Po 

In the similar fashion as in the univariate case they defined and <7^^ as 


Cpy]^ - '/77-i?7- 


F{N{U)) - 



1 

2 


l-T(jV(L)) | 

l-ai f 


( 6 . 21 ) 


for off-centered process and 


rM _ ( Pmn)) - Pmw F(N(T))-F(N(L)}\ 

^pTk — 


I 


1 


1 


/ 


2 ~ 0:2 2 ~ 

for off-target process respectively, where F(.) is the cumulative distribution function of X(.). 
They described the steps to be followed to calculate process capability index for a set of mul¬ 
tivariate quality characteristic data. 

1. Suppose Xi,X 2 , ■■■,Xn are n observation vectors come from a assumed multivariate distri¬ 
bution. 
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2. Transform the observations using N(X). Let Yi = N(X-,), Y? = N (Xo ) . Yn = N(X„). 

Here the distribution function F{y) = P(N(X) < y) is the structural function associated with 
X via Y = Nix'). Now NiX.) will possesses a distribution. 

3. Now assuming Yi,Y 2 , ...,Yn are observations, we will find a best fitted distribution for the 
observations using Q-Q plot technique or any goodness of fit test. 

4. Using this best fitted univariate distribution, PCI is to be calculated following the approach 
of Maiti et al. (2010). 

5. If there are more than one multivariate distributions assumed and corresponding structural 
functions are considered, then we will choose that structural function for which the transformed 
data set is best fitted. Here, to define generalized process capability index for multivariate data, 
conditional ordering has been used, and multivariate data has been transformed to univariate 
one using the concept of structural function. After transforming the data, it is as simple as 
applying generalized process capability indices for univariate data given by Maiti et al. (2010). 
It can be well understood and comfortably be used by the practitioners. 

7 Conclusions 

Process capability studies play an important role in the process control since they assist to 
decide whether a manufacturing process is suitable and the applications meet the necessary 
quality standards. The assessment of process capability, which now appears to be quite simple, 
involves significant dimensions in a practical set up. This is due to the fact that some of the 
conditions necessary to establish process capability not fully satisfiable. These conditions stip¬ 
ulate that the process has to be under control, that the process output has to follow the normal 
distribution and the observed values of the quality characteristics be statistically independent. 
Since these conditions are not fully satisfied in many manufacturing situations, the process 
capability analysis is becoming to be a critical issue. This situation was of course realized by 
many researchers and resulted in numerous publications in the literature. Therefore, process 
capability analysis is valid only when the process under investigation is free of any special or 
assignable causes (i.e., is in-control.). 

Recently, the research in the theory and practice of multivariate process capability indices has 
been very sparse in comparison to the research dealing with the univariate case. At present, 
for the multivariate capability indices, consistency in the methodology for evaluating this capa- 
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bility is not so much developed. Moreover, it is quite difficult to obtain the relevant statistical 
properties. Obviously, further investigations in this field are strongly desirable. In further 
correspondence, we will no doubt provide a fertile growth for new developments in the theory, 
methodology, statistical properties and applications of the PCIs. 
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